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The time-varying gravitational field produced by a Weber bar is used to explore math- 
ematical features of the linearized Einstein equation. We present a self-contained formal 
framework for the treatment of the linear field, which is applicable to several situations 



In this paper, the Weber bar is used as a specific example of a gravitational source, to illustrate 



o 
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^ where the standard quadrupolar formulas are not adequate. The expressions for retarded 

.^ integrals reveal a singularity associated with boundary conditions. Results from the trans- 

psj verse traceless subgauge are compared with the radiation calculated from retarded integrals. 

Lienard-Wiechert potentials are used in a treatment of the Weber bar as a collection of 

fj point particles and further possible applications are outlined. The Riemann tensor clarifies 

I the transition from near-field geodesic forces to tidal forces in the far field. 
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2^ I. INTRODUCTION 

o 

(N 

^^ various aspects of the mathematics and physics of the linearized Einstein equation. The discus- 

^^ sion comprises techniques for computing retarded integrals, including Lienard-Wiechert potentials. 

•^H These techniques are complementary to less detailed discussions of more general sources found 

X 

;-H in textbooks. For example, boundary conditions are a crucial aspect of the Weber bar, in that 

the outer few grams of material contribute as much to the gravitational disturbance as the inner 
hundreds of kilograms of mass. The Riemann tensor interpolates smoothly between gravitoelectric 
fields in the near field and tidal forces in the far radiation field. 

The aim of this work is to present a self-contained formal framework for the computation of the 
linearized field in situations which are different from those considered in standard textbooks. The 
typical approximation methods employed in the literature, based on the oscillating-quadrupole 
formula, involve several simplifying assumptions: far fields (radiation region), no retardation in 

* Email address: r31@psu.edu 

^ Email address: giovanni.modanese@unibz.it 



the source, symmetry of plane wave or spherical wave, TT gauge and absence of longitudinal 
field components. We relax all these assumptions. We also avoid the use of power series and 
Fourier expansions in the source, which might fail if sharp discontinuities are present. We provide 
an example of numerical and analytical computation of the near field. Within our framework, 
we describe the effects of the fields on detectors through standard formal ingredients (Riemann 
tensor, geodesic deviation, gravito-electric field components), but also considering possible forces 
and strains which are not typical of plane transverse waves. 

In Sect. [IT] the metric produced by a Weber bar is related to retarded integrals of the stress 
tensor. A mixture of Lagrange and Euler variables is used to establish notation. The technique 



for treating boundary conditions is shown explicitly. Section III describes the procedure used to 



express the quadrupole component of the radiation in the transverse traceless (TT) subgauge. 



Section III A discusses the (non) existence of a global gauge transformation between the harmonic 



gauge and the TT gauge. A local transformation can always be obtained, however. 



Section IV describes expressions used to calculate force gradients in a detector. Gradients in 
the gravitoelectric field and the second time derivative of the spatial components of the metric are 
treated as Newtonian forces. The Riemann tensor, which includes both geodesic and tidal forces. 



is introduced. In Sect.s IV A and IV B the response of a detector located on coordinate axes is 
expressed, resulting in the same response using either the retarded integral or transverse traceless 
metrics. 

In section |Vj expressions for near field components of the metric are evaluated, taking into 
account higher order terms in the retarded integral. Time-dependent metric components varying 
as 1/r^ and 1/r^ are computed. A numerical integration shows how metric components evolve 
from near-field to radiation field phenomena. 



In section VI the detector is located in a direction not on a coordinate axis. The descriptions 
of the forces in terms of the Riemann tensor are the same for the retarded integral and TT metrics, 
even though different components of the metrics contribute to Riemann tensor elements. The 
agreement between retarded integrals and the TT metric is important, since the derivation of this 
relationship in some texts ignores boundary conditions. 



In section VII, the Weber bar is treated as a collection of point particles, using Lienard-Wiechert 



potentials to express the gravitational field. The result agrees with that computed using density 



as a continuum distribution. In section VIII we present in a qualitative way some possible future 



applications of the method of Sect. VII, for instance to the computation of gravitational fields 



generated by stationary waves in a plasma. 



II. GRAVITATIONAL FIELD OF AN AXIAL OSCILLATOR 

The axial oscillator considered in this paper is similar to a Weber bar, which oscillates longitu- 
dinally in the z direction (Fig. 1). The linearized Einstein equation relates trace-reversed metric 
deviation elements and energy-momentum tensor elements as follows 
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Metric deviations are expressed in terms of the metric as follows 



9^J.uix) = ri^y + h^„{x); 7]^^ = diag (-1, 1, 1, 1) 



The trace-reversed metric deviation is defined as 
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h^uix) = h^y{x) - -ri^yh\{x) 
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and satisfies the harmonic gauge condition (analogous to the Lorentz gauge in electrodynamics) 
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The solution to the Einstein equation can be expressed as a retarded integral 
/i^-(xo,to) 
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For points outside the source, the fields (pi) can be evaluated using the series expansions 
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It is convenient to express the metric deviation as a sum of near and far fields. Time-dependent 
terms of first order in 1/ro correspond to radiation (energy transported to infinity), while terms of 
higher order refer to near fields. 

The energy-momentum tensor can be written in two diff'erent ways. In "Lagrangian" variables 
[2] one has 
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FIG. 1: Weber bar with longitudinal oscillations in the z direction. A generic point Q of the bar, with 
coordinate x, has at a given instant displacement w with respect to its equilibrium position, and velocity 
V. Both w and v are directed along z, and can be positive or negative. The retarded integral ([5]) gives 
the metric at an observation point P with coordinate xq . The upper and lower boundaries of the bar have 
coordinates z^op and ZBottom which vary in time. 



where p is the mass density, w a longitudinal strain variable and E the elastic modulus of the 
material. Alternatively, one can use "Euler" variables [3j 
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where p\ is a density perturbation (p = po + /^i)) ^ is the material velocity and Vs is the speed 

of sound in the material, given by Vg = 

continuity and stress-strain relations connecting -q^ , -q^ 

the energy-momentum tensor is equal to the mass current J^ multiplied by c, and the component 

T^^ is also called the pressure p. 

When the bar oscillates in its lowest frequency normal mode, the longitudinal displacement of 
any point of the bar with respect to its equilibrium position can be expressed as follows 



— . For an ideal elastic material, there are known 
and V at any point. The component T^^ of 
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where L is the length of the bar, e <C 1 is a perturbation parameter, and uj is the resonant frequency. 



Expressions for mass current and the pressure are straightforward. The pressure is written as 

p [z, t) = piVg = —E—— = —oj L pq— cos — - coswt (11) 

oz vr V L / 

The current density is expressed as follows 

dw T /"^^\ f-i^\ 

Jz [z,t) = p-^ = —epQiuL sm ( — I sm ut (12) 

In this section, only radiation is described, since the mathematics is simpler. Evaluation of near 
fields (which is impossible in the TT gauge) is postponed to Sect. W\ As usual in the calculation 
of the radiation field, the denominator |xo — x| in eq. (jsl) is replaced by tq = |xo| (distance from 
the center of the Weber bar to the observation point P). This is justified because the wavelength 
A = ^^ of the radiation is very long compared with the dimensions of the bar. In the following, 
we shall write explicitly this kind of approximation only when needed. 

The retarded integrals of the pressure and of the mass current are quite straightforward, because 
these quantities are of first order in e, so it is not necessary for their computation to take into 
account the movement of the upper and lower boundaries of the bar. That movement is itself of 
order e, and plays a role only in the computation of the integral of T^^ (see below). The retarded 
integral of the pressure is as follows 

h^^ixo,to) = ^ [ d^xp(x,to--) =-^ML'^u;^^cosiujto-kro) (13) 
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where M is the mass of the Weber bar. Evaluating the retarded integral of the mass current 
requires a position-dependent expression for the retarded time, 

/i°3(xo,to) = — ^ / d^3;j(x,to- ^°~""^ ) = J ML V n^ 4 cos (wto-^^-o) (14) 

roc^ J \ c / tqc* vr"^ 

where n = {nx,ny,nz) is a unit vector pointing towards the observation point. 
The expression for the retarded current contribution involves the integral 
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using the long wavelength approximation kL ^ 1. 

Evaluating the retarded integral of the mass density T^^ is more elaborate. The density is 
written as a sum of an initial constant density pQ and perturbed pi terms, where 
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The integral for h^^ has a hidden singularity, which is revealed by expressing only the time- 
varying component 
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where A is the cross-section area and the z coordinates of the top and bottom boundaries are given, 
as a retarded function of to) by 



(18) 



ZTop (to) = \L + eL cos (ujIq — kro + ^UzkL) 
ZBottom (to) = -\L - eLcos (ivto - kro - \nzkL) 

Here k is the wave number k = lo/c. 

The few grams of material in the boundary regions ZBottom{t') < z < —L/2 and L/2 < z < 
ZTop{t') contribute as much to the gravitational field as the hundreds of kilogram material in the 
interior. 

The time-dependent component of W^ consists of integrals over the interior and over the bound- 
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where $ denotes for brevity the phase $ = (wto — ^'^o) and kz is the z component of the wave 
vector k = kxi.. 

The boundary and interior terms have opposite signs, but comparable magnitudes. Their sum 
is 
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Note that the factor — ^ is of second order in the ratio of the length of the Weber bar (of the 
order of 1 meter) and the wavelength of gravitational radiation (of the order of 100 km). Also note 
that Weinberg ([3], p. 270) appends alternating positive and negative masses to the ends of the 
Weber bar. The mathematical justification for this procedure is unclear. The contributions to the 
metric from interior and boundary terms are comparable. 



Eq. (21) can be rewritten 
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In conclusion, the trace-reversed metric deviation from the retarded integral solution to the lin- 
earized Einstein equation in the harmonic gauge can be expressed in matrix form as follows 
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The non-trace-reversed form is 
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A metric of this kind is termed "Class I" by Boardman and Bergmann [3] . 
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III. TRANSVERSE TRACELESS METRIC 



In order to construct a spherical transverse-traceless metric we follow the recipe given by Finn 
[5]. This is based on the quadrupole moment Qij{t) of the mass distribution in the source {i,j = 
1,2,3) 

1 
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The tensor Qij (t) is combined with the projection operator Pjk (xq) = 6jk — rijUk, depending on 
the observation point xq, to obtain a tensor Qf-^ (xo,t) 



Qij (xo, t) = Pik{xo)QklPlj (xo) - -Pij (xo) QlmPlm (xq) 



(26) 



The transverse traceless metric hj^ (xq, to) is proportional to the second derivative of QJ^ (xq, t) 



with respect to t evaluated at the retarded time (to — ?^o/c) 
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For our oscillating Weber bar the tensor Qij{t) = Iij{t) — \lkk{'t)^ij has one time-dependent 
component, namely 
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Note that evaluating this integral requires the same boundary/interior treatment as that used for 
calculating h^^. 

After multiplying matrices and considering that the components /iqq and hf^ are zero due to 
the gauge condition, the result is 
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A. Gauge transformation 



There exists no rigorous proof that the TT metric (29) can be related to the harmonic metric 



(23) through a gauge transformation valid in all space. It is well known that the TT gauge can 



be imposed on a plane gravitational wave, leaving only the physical transverse components of the 
field. The same can be done for fields which are superpositions of plane waves, but the TT gauge 
is not obtainable in general, as we show in the following. 

The standard argument proving that it is possible to pass from the harmonic gauge to the TT 
gauge outside a source goes as follows [I\. Suppose to have a trace-reversed metric h'^'^{x) which 



satisfies the harmonic gauge condition d^h 



flU\ 



and the field equation d^'dah^^ix) = 0. The 



infinitesimal gauge transformations of the coordinates have the form 



x'f" = x'' + Cf'{x) 



Under these transformations the metric changes as 



Kui^') = V(^) -^f^uix) 



(30) 



(31) 



where £,^u[x) is defined as 



iiiv{x) = d^iu(x) + du^^{x) - Vfiudp^^ix) 



(32) 



The four-divergence of the metric changes as 

(a'^ V(:e))' = d^h^,{x) - d'^do.U^) (33) 

therefore the harmonic gauge is not spoiled in the transformation, provided d'^daCuix) = 0. By 
suitably choosing four functions S,u{x) which satisfy this latter condition, and inserting them into 



(31), (32), we can hope to cancel the trace and some components of /i^'^(x), eventually obtaining 



the TT gauge. Relation (31), however, is not just an algebraic equation, but a partial differential 
equation (except in some simple cases, like for plane waves). This introduces some complications. 
Usually one starts by requiring that the trace of h '^'^(x) is zero and then attempts to change 
the h ^"^ components. The null-trace condition translates into the following condition on £,u{x) 

r]^%Ax) = => r/''" V(x) - v^''C^.u{x) = ^ 2dpeix) = -n^'h^Ax) (34) 

So, in order for h^'^(x) to be traceless, the gauge transformation must satisfy the equations (in 
units such that c = 1) 

' do'C,{x) - Ae^(x) = 
5oe°(^) - d^eix) = -h^'hf.Ax) 
The solutions of these partial differential equations must respect some conditions at the boundary 
of the source, because inside the source it is impossible to pass to the TT gauge {d^daS,u{x) = is 
not true there. See Fig. [21) 



Even disregarding the boundary conditions, it can be impossible to find a solution of (35) in a 
general case. In the wave region the trace of the (given) harmonic metric has time dependence of 
the form ri^^"h^u[x) = /(x)e*'^*, so we make the ansatz ^p(x, t) = Xp(x)e*'^* and the second equation 



of (35) becomes an equation for Xp{^) 



iwXo(x)-aaXx) = -^/(x) (36) 

We can solve this by requiring that 9jX*(^) — ^^'^ Xo(x) = ~2iu)fi^)- There is a trivial solution 
X*(x) = , but this certainly does not allow to cancel the components /iQi in the next step. If we 
take x*(x) to be a transverse plane wave with wave vector k we can satisfy the condition 5iX*(x) = 



and then also the second equation of (35), provided |k| = uj; but after that, we will be able to 
cancel the components liQi only if they are plane waves themselves. In spherical symmetry we can 
satisfy the condition diX^{x.) = with functions x*(^) 7^ which are like the components of an 
electric field in vacuum and behave like 1/x^; but at the same time these functions should satisfy 



the second condition of (35), namely — a;^x«(x) — d^dkXii'^) = 0) and this is clearly impossible. 
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FIG. 2: Local implementation of the TT gauge. At any point P in the radiation region the gravitational wave 
can be locally approximated with a plane wave of vector k. The components Xi of the gauge transformation 



to the TT gauge (see eq. (36l) are then chosen of the form Xi(x) = 9ie , with 9iki = 0. In order to obtain 



a global implementation of the gauge condition, however, we should solve eq. (36) over all space outside the 
source S, with suitable boundary conditions at the border B of S. 

Nevertheless, the gauge transformation holds locally because it is always possible to approximate 
the wave with a plane wave in a region much smaller than its wavelength. This provides an indirect 
proof that the Riemann tensors for the harmonic gauge and the TT gauge are the same at any 



point (an explicit calculation is given in Sect. VI). This ensures, in turn, that the physical effects of 



the field on detectors are the same, as expected. Nevertheless, in the following sections we analyze 
in more detail the behavior of gravitational waves detectors also in more general situations, where 
the TT gauge is not applicable. 

IV. GRAVITATIONAL FORCES ON A DETECTOR: BASIC DEFINITIONS 

Gravitational forces on a detector come in two varieties. For observation points in the near field, 
the metric varies over distances comparable to the distance from source to observer. In the far field, 
the metric varies over distances comparable to the wavelength of radiation. For measurements in 
the near field, gravitational forces can be expressed in terms of the geodesic equation 

d^x^' ^„ dx^dxP 

^'u,-j--j- (37) 



dT^ ""P dT dr 
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where x^{t) describes the trajectory of a test particle, and the Christoffel symbol describing vari- 
ations in the metric are defined as 

Ku = 5'"'' ki^, p] ; [criy, p\ = 2 ^3up,u + gpa,u - gau,p) (38) 

Forces on a stationary test particle can be described in terms of a gravitoelectric field 

El = -r,,c' = ^^ + ^T^ (39) 



At larger distances from the source, tidal forces associated with time derivatives of the spatial com- 
ponents /i*-^ are important. The geodesic distance between two points xq and x in a detector differs 
from the coordinate distance. Since electromagnetic forces and wavefunctions are constrained to 
propagate along the geodesic distance, the two points must have a relative acceleration to balance 
stresses and strains 

^ = i^X^- (40) 

where X = xq — x is the coordinate separation of the two points. The transition from near field to 
far field behavior is expressed smoothly in terms of the geodesic deviation equation, which involves 
the Riemann tensor, 

where D denotes a covariant derivative, 6x denotes the coordinate difference between a spacetime 
point x^ and a reference point. A convenient formula for the Riemann tensor in covariant form is 
given by 

Rpapiy = [crj^, p],^ - [o-^i, p],^ + Kp [P^> tt] - Tl^ [pfi, vr] (42) 

Ignoring the higher order terms in the Riemann tensor, the relative acceleration between neigh- 
boring points can be approximated as follows 

D^6x^ 



r\j 



R'q^qC^Sx^ (43) 



D^2 - ^^OjO 

The four-velocity of a point is approximated in terms of the time-like component 

dx^ 

-^ = 'yic,Vcc,Vy,v^)^{c,0,0,0) (44) 

The velocity of the detector is assumed to be small compared to the speed of light. 
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In this work the relation between gravitoelectric and tidal forces is displayed explicitly, for the 
radiation field of the Weber bar. In the linearized approximation, using the expression for the TT 



metric, the gravitoelectric field defined in (39) is given by the vector 
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2 ^ 2 
It is not obvious at this stage that the Riemann tensor nullifies the gravitoelectric field in radiation. 



Gravitational wave detector on the x axis 



Having explicit expressions for the harmonic metric ( 23 ) and the TT metric ( 29 ) allows gravi 



toelectric and strain field forces on a detector to be evaluated. The detector consists of two rings 
centered on the x axis. The rings each have a radius 5, oriented parallel to the yz plane. In the 
absence of gravitational radiation, the distance between the rings is s. Anticipating the need to 
evaluate phase gradients, define the origin of coordinates as the center of Ring 1 (see Fig. pi). For 
an observer on the x axis, the TT strain field is given by 
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/ij • {ux = 1) = —H cos {ujt — kri — kx) 
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Consider points in rings 1 and 2 originally at the locations 

xoi = (0, B cos (j), B sin cj)) 
xo2 = {s, B cos cj}, B sin cf)) 



The tidal forces (eq. (40)) produce the following accelerations 

/o o\ / ^ 
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(50) 
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FIG. 3: General detector for gravitational waves. The detector consists of two rings centered on the x axis. 
The rings each have a radius B, oriented parallel to the yz plane. In the absence of gravitational radiation, 
the distance between the rings is s. 



In particular, the longitudinal force is zero. For the harmonic metric, the strain field is 

^ -1 0^ 

h^j {ux = 1) = —H cos {ujt — kri — kx) 
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The tidal forces produce the following accelerations 
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The tidal force includes a longitudinal component, absent from the TT force. In addition, the 



harmonic metric includes a gravitoelectric field from eq. ( 45 ) 

1 
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The phase difference causes a relative acceleration between the rings. We denote this acceleration 
with the superscript "(/" , in order to distinguish it from the "strain" accelerations computed above. 
The superscript ".ff" means, as above, that this acceleration is computed in the harmonic gauge. 

X2 '^ — x^ '^ = —^Hujc (sin {ujt — kri — ks) — sin {uit — kri)) 

r / 'Al (^^) 

= \H up' s cos {oot — kri — ks) — \HijJC \{ks) sin {oot — kri) + O [{ksY 



The relative acceleration produced by the gravitoelectric field (55) cancels the longitudinal com- 



ponent of the tidal force (eqn. (53)), up to second order in the phase ks. 

Both the harmonic and TT metrics produce transverse motion of the detector rings and sur- 
rounding instrumentation. The harmonic metric also produces a longitudinal oscillation of detector 
and instrumentation. However the sum of the longitudinal forces from the gravitoelectric field gra- 
dient and strain fields cancel in first order in the phase ks, producing a small net oscillation in the 
separation of particles. The expression of T^^ in eq.s ([8|, (19]) for stress in a material requires a 
strain gradient. Since only relative motion is detectable, the harmonic metric produces almost the 
same detector response as the TT metric. 

The geodesic deviation equation relates strains in the x direction to Riemann tensor elements 
of the form Riqjq. 

^1010 = [00, 1] ;l - [01, 1] = --ffoo,ii - ^^ii.oo (56) 

[00, 1] = 2 (^oi'O + 510,0 - 500,1) = -2500,1 (57) 

[01, 1] = 2 *^^i^'° + ^i°'i ~ ^°^'i) " 2^^^'° ^^^^ 

^1020 = 0; i?io30 = (59) 

For the harmonic metric 

g'oo = —511 = —H cos {ut — kri — kx) (60) 

Hence 

(p-n (p-n f uP'\ 

500,11 + 511,00 = "o-^ + ^^ = ^cos (wt - kri - kx) [ -k^ + — j = (61) 

since w = kc. 

The derivative 500,1 is the x component of the gravitoelectric field. The gradient of the grav- 
itoelectric field cancels the second time derivative of /in. Since the Riojo terms of the Riemann 
tensor vanish, the x component of the strain vanishes. 
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B. Gravitational wave detector on the z axis 



Consider the detector in Fig. 3 located on the z axis. At this location, the TT metric (29) 



vanishes, and the harmonic gravitoelectric field (45) is non-zero. A longitudinal component of the 



harmonic tidal strain field persists, resulting in oscillations between the two rings. 



-H, strain 



Z 



■H, strain 
1 



-Hlo' 



f 0^ 


/ 







Voo ij 


V 



i?coS( 
Ssinc 



\ 



/ 



cos {ojt — kri — ks) 



(62) 



The effect of the gravitoelectric field on the difference Z2 — z\ is as follows 

^2 '^ — z^'^ = ^Hojc (sin {ojt — kri — ks) — sin {ujt — kri)) = 

r / -Al (63) 

— \Huj'^s cos {ojt — kri — ks) + \Hujc \{ks) sin (wt — kri) + O ( {ks 

The gravitoelectric field produces longitudinal motion of the detector and surrounding instrumen- 
tation, which cannot be detected with local measurements. The gradient of the gravitoelectric field 
cancels the tidal strain, up to second order in the phase ks. The relevant Riemann tensor element 
is 



-^3030 — 530,03 — ^500,33 



1 



:533,oo — 



(64) 



since 



500 



-530 = 533 = -H cos {ujt - kri - kz) 



(65) 



The combination —930,0 + 5500,3 is the gravitoelectric field. The gradient of the gravitoelectric field 
cancels the second time derivative of 533. Since the -R3030 Riemann tensor element vanishes, the 
strain force in the z direction is zero. 

V. TIME-DEPENDENT NEAR FIELDS 



One technique for computing near fields from retarded integrals ([S]) is to use series expansions. 
The denominator can be expanded as in eq. (p]). The time retardation involves the distance between 
source and observer, eq. M. The result is a double infinite series of terms, which can be evaluated 
using trigonometric identities and integration by parts. Retaining only time-dependent terms up to 
second order in wavenumber k and third order in l/rg, the trace-reversed metric deviation tensor 
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is as follows 



/i^'^-(xo,to) 



SGML^e 



TT^C^ 



_fe2 (l-3n2) 






\ 



— ^ cos $ + % sin ^ 



ro 







V =^ cos $ + % sin $ ^ cos $ 

where # = (wto — kro). As an example, the 03 metric element is related to the mass current 



(66) 



/ 



hnear (^Oj *o) — ^"^ / d X J ( X, to 



r^r^ 



n • X 



The equation involves the following integral 

L/2 



t03 



dzsin f -— ] n-yZ 



\ L 



L' 



vr^ 



^Uz [svciu — ucosu 



7r/2 
-7r/2 



2L' 



-n. 



vr^ 



-L/2 



A. Gravitoelectric and gravitomagnetic fields 



(67) 



(68) 



Forces on a test mass can be expressed in terms of the geodesic equation 



dfi 



Hnc^ - rUv^ - TLv'' 



oo<- ~ -^ Oj 



kO^ 



Remembering that 



1 



dgoi _ dgoj\ 
2 ^ V dx^ ' dx^ dx^ J 



n, = ^9'' i ^ + 



J^ 00 — off 1 ^^^ + 



2 \ Sx*^ ^x" dx^ 



(69) 



(70) 



normalization and sign conventions are chosen so that the force on a test mass can be expressed 
as follows 



F = mtj„ + mw X Bg m ^ v Sj 

2 ot 



(71) 



where ij is a unit vector in coordinate direction j. The gravitoelectric field is expressed in terms 
of metric elements as follows 



df^2^ Ikf ~ A^ Ikf ^ 4^ Ikf ^ ^^df''^ 



El = -c'no = c^^ + ^c 



+ 4' dx^ ^"^ dt '^3 
It is convenient to express the gravitomagnetic field as the curl of a vector potential 



(72) 



Bg = V X h = c (-ay/^°^ 9^/l°^ 0) 



(73) 
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To conform to the geodesic equation, the vector potential must be chosen in terms of covariant 
metric elements 



h = C (/loi, ^02, M = -C (/l°\ /i°', /i°') 



(74) 



For bookkeeping purposes, derivatives are labeled according to orders in k and 1/r. The most 
important terms in the near field region involve gradients of the l/rp term of the h?^ metric 
element 

2GML'^e r(l-3n2) 



^9 ( ^4 ) " 7r2c2 



COS <I>V 



(75) 



E 



9 \ ^4 



-y- COS <l> [Axz^ - y^x - x^, Ayz^ - x^y - y^, 2z^ - 3x'^z - 3y'^z) (76) 



Ti'^c^r' 



The leading term for the gravitomagnetic field involves the curl of the spatial factor in the 1/rQ 



term of the h?^ metric element 



B 



9 \ ^3 



SGML^e 



A;sin<I> 



V X 



zz 



24.GML^eksm^ ^ 

Y^ ^i-yz,xz,0) 77) 



B. Displacement current 



The displacement current involves a relation between the curl of the gravitomagnetic field and 
time derivative of the gravitoelectric field. The leading term (linear in k) in the curl is 

V X B„ ( -r ) = n o V ^ (x^ + y^x - Axz'^,y^ + x^y - Ayz^ , Zx^z + 3y^z - 2z^\ sin $ (78) 



For comparison, the time derivative of the gravitoelectric field ( 76 ) is 



^ ' ' =^ ' ^sin$(4rEz2-y2j;-x^,4yz^-x^y-y^,2z^-3x^z-3y^z) (79) 



vr^c^r' 



ct \r^ 
The two expressions are related: 



„ ^ , k\ ,dE„ f k 



(80) 



The relative signs in the displacement current equation are the same as the corresponding equation 



in electromagnetism. However the relations between fields and sources (eqns. (72) and (73)) have 
opposite signs to the corresponding relations in electromagnetism. 
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C. Faraday induction 

Faraday induction involves a relation between the curl of the gravitoelectric and time derivative 



of the gravitomagnetic field. Since the curl of the gradients in eq. (72) vanish, only the time 
derivative of the vector potential contributes to the curl of the gravitoelectric field. 

Eg = V (term) + c— — (81) 

The equation for the gravitomagnetic field is 

Bg = V X h (82) 

Since space and time derivatives commute, 

VxE, = -^ (83) 



Note that hz = cho^ = —ch^^, resulting in the minus sign in eq. (83). Conclusion: both Faraday 
induction and the displacement current obey laws similar to those for electromagnetism, in the 
near field region. In vacuo the equations have the same sign structure as in electromagnetism. 
However the relations between fields and sources have opposite signs for gravitoelectromagnetism. 



What about the h^^ derivative forces? The force equation ( 71 ) includes velocity-dependent 
forces on a test mass associated with time derivatives of spatial components of the metric. The 
three diagonal elements of the metric deviation are as follows 

/jll = /j22 = 1 (/,00 _ ;^33) . ^33 ^ 1 (/^OO ^ ^33) (34) 



All three elements include the term 



SGML'^e (1 - 3nl) 



h''= ^ , „ ^ ^cos$ + ... (85) 



The time derivative of these elements is of order k/r^ 

dh'' / k\ SGML'^e [l - 2,nl) 



k sin $ (86) 



cdt \r^ J 2Tr^c'^ r'^ 

These terms are of the same order (linear in k) as that of the leading term ( 77 ) in the gravitomag- 



netic field. 

Concerning the relation between the harmonic gauge and the TT subgauge, remember that the 
linearized Einstein equation (nj) is derived by applying the harmonic gauge condition d^h'^'^ = 0, 



19 



to disentangle relations among metric elements. As a consequence, derivatives of metric elements 



(66) should be related; for example 



8W 8ft»3 

The z derivative of the spatial factor in the l/rg term of hP^ results in a l/rg term. This term 
is matched by the time derivative of the phase factor in the 1/rg term of h?^. The Lorenz gauge 



condition requires relations among terms with differing powers of l/rQ. The radiation metrics (23) 



and (29) lack the higher order terms to fulfill the harmonic gauge condition. These metrics can be 
considered as solutions to the Einstein equation only in a restricted sense. The metrics must be 
treated as if they were locally plane waves, so that derivatives of the 1/ro factor are to be ignored. 



The mathematics of this ansatz is described more carefully in section III A 



D. Numerical integration 

Another technique for evaluating retarded integrals is numerical integration. As an example, 
the behavior of the gravitational field along the z axis was computed using numerical integration. 
Specifically, the integrals calculated numerically are as follows 

hZm{z,M= I dz^P^=Y,Az^^^^^^^^^ (88) 

J po — z "!— • Zo — Zi 

^Bottom 

No series expansions were used, so that the main source of error is numerical precision, rather than 
truncation of a series. 

The retarded time t' assigned to each z value in a coordinate grid is t = to — (zq — z)/c, with 
to chosen as one tenth of an oscillation period. The Weber bar was divided into 200 grid cells. 
The perturbation parameter was chosen as e = 0.005. For sufficient accuracy in computing h^^^ 
(interior), a Simpson rule integration over the width of a cell was used. In Visual Basic code. 



H00=H00+H00cof*(RHiiiii+4*RHcn+RHpl)/6, where RHcn is the mass density perturbation (16) at 

the center of a cell, and RHmn and RHpl are density perturbations at the bottom and top of a cell. 

Showing the transition from near-field to far-field behavior for the real Weber bar is problematic. 



because the cancellation between h^^ (interior) ( 19 ) and h^^ (boundary) ([2|) requires high numerical 



precision. To obtain a numerically significant result, the frequency was arbitrarily increased by a 
factor of one hundred over the 1667 Hz resonance frequency of the real Weber bar (see Fig. HI). A 
coefficient of GML'^ /c'^ was set to unity. 
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0OE*CO -■- h03 



FIG. 4: Metric components in the transition region between near and far fields. At a distance of one wave- 
length, h^^ merges smoothly with Jv^^ and h^^. Close to an oscillating source, the time- varying gravitational 
disturbance is measurable. For example, a torsion pendulum is used to search for deviations from Newtonian 
gravity 6J. Beyond a distance of a wavelength, ft,'^°, hP^ and h^^ oscillate in unison, resulting in zero for the 
-^3030 Riemann tensor element. 

VI. RIEMANN TENSOR FOR ONE WEBER BAR IN OFF-AXIS DIRECTION 

Interesting Riemann tensor elements are required, for a direction not on a coordinate axis. 
Forces on a detector located in the direction 



n= (0.6,0,0.8) 



(89) 



are evaluated using both the retarded integral and transverse traceless radiation metrics. For this 
direction, the harmonic and TT metrics are 

/o ^ 

0.1152 -0.0864 



h^""'^ = H 



/^0.82 0.8 ^ 

-0.18 

-0.18 

\ 0.8 0.82 y 



h'"'- 



TT 



H 



-0.18 
\0 -0.0864 0.0648 / 



where 



H 



8GM 



L eoj cos {ojt — kr) 



(90) 



(91) 
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Harmonic TT 



R-ioio - Iffoo, 11 - Isii.oo = -0.0576 -|5ii,oo = -0.0576 

hO.5 • 0.82 • 0.62 _,_ Q 5 . Qjg] [_Q 5 . 0.1152] 

R1030 -5500,13 + 5503,10 = 0.0432 -^531,00 = 0.0432 

R.2020 -5522,00 = 0.09 -5522,00 = 0.09 

R.3030 530,03 — 5500,33 — 5533,00 = 0.0324 —5533,00 = 0.0324 

R.3131 —5511,33 — 5533,11 = —0.09 513,13 — 5511,33 — 5533,11 = —0.09 

R-3130 5.903,13 — 5533,10 = 0.054 5531,03 — 5533,10 = 0.054 



TABLE I: Elements of the Riemann tensor computed in the direction of eq. (89). A common factor oioj^H" 
is omitted. The numerical values of the single terms are given explicitly in square brackets only for i?ioiOj 
as a calculation example. 

Below is a table of Riemann tensor elements. In order to evaluate the second derivatives, express 
the coefficient H as follows 

H = 7: — jL^euP' cos iojt — kri — kn ■ x) (92) 

Then 

H^o = ujH'; H^i = -con^H'; H^s = -con.H' (93) 

H,oo = io^H"; H^oi = -uj\,H"- H,o3 = -u;\,H" (94) 

H,u=io\lH"; H,ri = io^n,n,H"; H^^s = u^nlH" (95) 

Augment the table of Riemann tensor elements with numbers, omitting a common factor of uj'^H" . 
Different components of the harmonic and TT metrics appear in the Riemann tensor. However 
the combinations "magically" end up with the same Riemann tensor elements (Table ^ . 



Incidentally, the gauge transformation (30) relating the two metrics in eq. (|90[) is 



U''^] = [x^] - ^-^L^ew (0.41, 0.246, 0, -0.472) sin [ujt - kr) (96) 



TT^rc^ 



VII. TIME- VARYING FIELDS FROM POINT MASSES 

An oscillating Weber bar can be considered as a collection of moving point masses. A sum- 
mation over point mass pairs produces metric elements consistent with results using a continuous 
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distribution of masses, provided a correct interpretation of the retarded Greens function is used. 
Consider point masses located at the top and bottom surfaces of the Weber bar. The locations of 
the masses as a function of time are as follows 

XTop(i) = (0,0, -L + eL cos (wt)j (97) 

^Bottom (t) = (o,0,--L-£L cos (ujt) j (98) 

The configuration can be described in terms of a conserved stress tensor, 

TOO (x, t) = mc^6 (x) 5 (y) [5 {z - ztop (t)) + S {z - ZBottorn m (99) 

r°3 (x, t) = -mc5 (x) 6 (y) [5 (z - ZTop (t)) - 5 (z - ZBottom {t))] eLuj sin wt (100) 

r33 (x, t) = mLu^ed (x) 5 (y) 9 {zTop (t) - z) 6 {z - ZBottom it)) cos wi (101) 

The T^^ element can be interpreted as a massless spring, with a tension Tspring = mLuj'^e cos ujt. 
The W^ metric element must be expressed as follows 

/,oo (^^^ to) = ^ [ dt' I rf3^ ^°°(^'^P ^ (t^ _t'_ ^^i^^ (102) 

c^ J J |xo - x| V c J 

The standard technique for computing Lienard-Wiechert potentials requires to perform first the 

integration over d^x. Briefly, the dt' integration variable is changed to the resulting argument of 

the delta function 

c 
The change in integration variable results in a velocity-dependent factor multiplying the retarded 
value of the integrand 



di^ _ 1 - 1 + ^Top • "VTop f rxop ■ ^TTop 'Y I'v^ 

(It I _ fTop-VTop c V c y \c^ 

c 

where 



(104) 



xq - :>^Top {t d^Top {t . . 

^Top = -, 77777; ^Top = 77 (105 

|xo-XTop(r)l dt 

Similar expressions apply to the point masses on the bottom. Summing the two contributions 

AG 

/l°° (xo, to) = 4 [2mC^ + mc {fxop ■ ^Top + r Bottom ■ ^Bottom)] ^gj (106) 
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An iterative procedure is used to compute the retarded times, which are sUghtly different for the 
two masses. As a first approximation, 



|xo| rQ 

tret = to = ^0 

c c 

The distances of the two masses from the observation point at the retarded time tr are 



(107) 



rTop = |xo -XTop{tr)\ 



Xn Lz 

" 2 



ro - -Lriz 



(108) 



fBottor, 



xo + -Lz 



ro + -Lu;, 



(109) 



The velocities of the two masses at the retarded times are computed using an iteration of the 
retardation times 



,/ ' lop rQ Luz I 

tTop — ''0 —to 1 — z — , I ; 



rTop _ rQ LUz , _ r Bottom _ . ro Ln^ 

— to H -T, ; t^ ^^ — to —to - 

c c 2c c c 2c 



(110) 



so that 



^Top = -sLuj sin ( ujto - kro + -k^L \ z; ^Bottom = eLio sin ( uto - kro - -k^L ] z (111) 



For brevity, denote the phase at the origin as $ = coto — kro. The unit vector rTop in eq. (106) is 
as follows 



Xo — ^Top ^ Xo - ^Top I -, , XTop • n 

rTop = ■; r = I 1 + 



l + ^^?:^^U-^^ = n + 0(^) (112) 



|xo - XTopI ro \ ro J \ ro J tq 

and similarly for f Bottom- The distinctions among rTop, rBottom and n are ignored. 



The time dependent component of h^^ in eq. (106) involves the sum 



Vsum = n • (vTop + ^Bottom) = n^eujl (- sin ($ + \kzL) + sin ($ - \kzL)) 
= -2nzeLoj cos ^ sin [^kzL) = -n^e^L^cos^ 



(113) 



Hence 



4G 



/i°° (xo,to) = — 4 [2mc^ - mw^L^n^e cos $] 



roc 



(114) 



The retarded integrals of T and T can now be computed, ignoring the velocity-dependent factor 



ofeq. (104) 



/i^^(xo,to) 



AG 

roC^ 



-mceLu) 



sin ( $ + -kzL \ - sin ( $ - -k^L 



AG 

roc^ 



mceL CO n^ cos ^ 
(115) 
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Hence 



AG 



h^''{^o,to) 



h^^ (xo,to) = J (mceLuj'^) Lcos$ 

roc** 

f nl n^\ 





V n, 1 / 

A. Sum over mass pairs 



(116) 



AG 



2 t2 

rribj eL 



cos {(jjtQ — kro) 



(117) 



The general relation between the location Zi of a point mass and its original position Zj is 



Zi = Zi + w = Zi + eL sin ( -^ ] cos u)t 



irzi 



L 



Sum (integrate) Z^ over all mass pairs, 

"ITZ. 



Zf = { Zi + eLsm(^] cos ujt^" ^ ~^ 



'vrz,- 



z.i + 2zi£L sin ( — ^ ) cos wt 



The time dependent part of the second moment is 

Qzz (t) = y 2miZi£L sin I — -^ 1 cos cot 



pairs 



Replace nii by rrii = MAzi/L, to convert the sum (120) to an integral 

L/2 



Qzz (t) 



M 

T 



^ TT 2' \ T 

dz2zeL sin ( —— ] cos cot = AeM—^r cos cot 



L 



vr^ 



-L/2 



hf"" (x, t) 



G 



The sum over mass pairs gives the result 

/ n^ n^ \ 





V n, 1 / 

The result agrees with that computed using density as a continuous distribution. 



vr^ rc'^ 



ML^eoj"^ 



cos {ujt — kr) 



(118) 



(119) 



(120) 



(121) 



(122) 



VIII. FURTHER POSSIBLE APPLICATIONS 



As we mentioned in the Introduction, the computational scheme presented in this paper can be 
applied to physical sources for which the usual quadrupolar formula is not adequate. 
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Examples of sources which do not have a symmetry compatible with the TT gauge are the arrays 
of mechanical or electromagnetic oscillators devised for the laboratory generation of high-frequency 
gravitational waves (HFGWs). Such generators have been discussed, among others, by Grishchuk 
[7]. They are still far from concrete realization, but it has been estimated that they could generate 
HFGWs with amplitudes only 100 times smaller than the main astrophysical sources. 

Retardation effects in the source occur when the wavelength of the gravitational wave is com- 
parable to the size of the source, so that the waves generated at different points do not add in 
phase. A typical example are electromagnetic sources made of resonant cavities (see [8j for an exact 
solution with a thoroidal cavity). One can also conceive an intermediate case, i.e. generators which 
contain stationary waves of a material medium, having propagation velocity much larger than the 
sound velocity of compression waves in Weber bars. Also in this case the retardation effects are 
important. 

Consider for instance a plasma supporting stationary waves of the Alfven type. These waves 
can have propagation velocity close to c and frequency up to several MHz. Their wavelength is 
then such to produce retardation effects in the source. It is interesting to observe that in this 
case the emitted gravitational wave can have an oscillating-dipole component without violating 
the conservation of momentum in the source. 

Namely, consider in the stationary plasma wave of wavelength A5 (red in Fig. pi) two points A 
and B, which are opposite with respect to the center of mass of the plasma. The plasma particles 
at these points move in opposite directions (red arrows), thus keeping the center of mass at rest. 
Each of these elementary oscillators generates dipolar waves with opposite phases (black and blue 
lines in Fig. l5| and wavelength Ac- Such waves would interfere destructively if their wavelength 
was much larger than the distance between A and B] but in this case the interference is not entirely 
destructive. 

Note that the "dipolar component" generated in this way has an amplitude comparable to the 
quadrupolar component, because in general the amplitude ratio of the two components is of the 
order of d/Xc where d is the size of the source. In fact, the multipolar expansion should be 
replaced by the computation of a full retarded integral. 

IX. CONCLUSIONS 

This document uses "pedestrian" but comprenhensive techniques to show a connection between 
the Lorentz-harmonic (retarded integral) and TT (transverse traceless) metrics. The local and 
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FIG. 5: Elementary illustration of retardation effects in a source of gravitational waves (plasma with high- 
speed transversal stationary waves). Two mass elements A and B of the plasma oscillate in phase opposition 
(red), so that their center of mass is at rest. The dipolar waves generated by A and B (blue and black) 
would cancel if their wavelength \q was much larger than the length A5 of the stationary "matter wave" (as 
it usually happens). In this case, however, the two wavelengths are comparable and the negative interference 
(here only shown in the z direction) is only partial. Note that in reality the near field is different from a 
pure sinusoidal function; the figure only shows that the argument for the total negative interference of the 
dipolar components is not applicable. 

global aspects of this connection are emphasized. The harmonic metric includes a propagating 
gravitoelectric field, which is absent in the TT metric. Gradients in the gravitoelectric field com- 
bine with tidal forces from velocity dependent terms in the geodesic equation, resulting in purely 
transverse local disturbances in the far field. The concepts of gravitoelectric and gravitomagnetic 
fields are valid in regions wherein time derivatives of spatial components of the metric are negligi- 
ble. A cross-over from near-field geodesic forces and far-field tidal forces occurs at a distance of a 
wavelength from a time- varying source. 

The physical effects of the metrics are evaluated for a general kind of detector composed of two 
parallel rings, which are displaced and deformed by the field. The harmonic metric is computed 
also in the near-field region, both analytically and numerically. The analogues of the displacement 
current and Faraday induction are displayed. A suitable application of the retarded Green function 
allows to re-obtain the harmonic metric of the oscillating Weber bar as an integral over point masses. 
Some possible further applications of this technique are discussed, e.g. to the field generated by 
a stationary wave in a plasma, which is expected to exhibit relevant effects of retardation in the 
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source. 
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